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ABSTRACT 

We present structural parameters for 204 stellar clusters in the Small Magellanic Cloud derived from 
fitting King and Elson, Fall, & Freeman model profiles to the V-band surface brightness profiles as 
measured from the Magellanic Clouds Photometric Survey images. Both King and EFF profiles are 
satisfactory fits to the majority of the profiles although King profiles are generally slightly superior 
to the softened power-law profiles of Elson, Fall, and Freeman and provide statistically acceptable 
fits to ^90% of the sample. We find no correlation between the preferred model and cluster age. 
The only systematic deviation in the surface brightness profiles that we identify is a lack of a central 
concentration in a subsample of clusters, which we designate as "ring" clusters. In agreement with 
previous studies, we find that the clusters in the SMC are significantly more elliptical than those in 
the Milky Way. However, given the mean age difference and the rapid destruction of these systems, 
the comparison between SMC and MW should not directly be interpreted as either a difference in 
the initial cluster properties or their subsequent evolution. We find that cluster ellipticity correlates 
with cluster mass more strongly than with cluster age. We identify several other correlations (central 
surface brightness vs. local background density, core radius vs. tidal force, size vs. distance) that can 
be used to constrain models of cluster evolution in the SMC. 

Subject headings: globular clusters: general — galaxies: star clusters — galaxies: Small Magellanic 
Cloud 



1. INTRODUCTION 

The study of the structure and evolution of stellar clus- 
ters has been reinvigorated by the dis covery of young, 
star clusters in ongoing gala xy mergers l|Holtzman et alJ 
119921 iWhitmore et al.1 119931) and poststarburst systems 



ijYang et alJ l2004|) . The desire to understand the for- 
mation and evolution of these objects, and to use these 
objects to provide a measure of the time and efficiency 
of star formation during mergers, provides yet another 
connection between studies of distant and Local Group 
galaxies. The Magellanic Clouds contain numerous 
young clusters that can be studied in detail because of 
their proximity and can guide our understanding of clus- 
ter evolution. A first step in the study of the cluster 
population is the measurement of such physical charac- 
teristics as the surface brightness profiles. We present 
surface brightness profiles and model fits for 204 clusters 
in the Small Magellanic Cloud. 

Our understanding of the surface brightness profiles 
of star clusters is dominated by two particular mod- 
els. The first generically successful set of models was 
empirically derived to reproduce the surface bright- 
ness profiles of Galactic globular clusters (|Kind [19621. 
These have become known as King profiles and were 
later derived from tidal ly limited models of isother- 
mal spheres ijKinel H"96fifl . Consequently they are ex- 
pected to accurately describe bound stellar systems, 
if they have isothermal and isotropic stellar distribu- 
tion functions, stars of the same mass, and reside 
within a tidal field exerted by a nother objec t. Ex- 
tensions to non-isotropic systems (|Michidll963j) . rotat- 
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ing systems l|Prendergast k, Tomeil 1197(1 IWilsorJ | 1975[) . 
anisotropic and rotating s ystems llLupt on fc Gunr J l987). 
and multimass model s l)Da Costa and Freeman! |1976t 
iGunn and Gri ffin 1979) exist. Only clusters in the near- 
est galaxies, including the Milky Way of course, are suf- 
ficiently well-resolved to test whether King models, and 
the corresponding extensions, are proper descriptions. In 
general, the basic King model is sufficient to explain the 
observed profiles of Galactic globular clusters, typically a 
surface density or surface brightness profile, although the 
extensions have been tested in a few clusters with more 
detailed obs ervations (for examp le, using radial veloci- 
ties in M 13 l)LuT)ton et alJ ll9871). As such, they can be 
viewed either as highly successful fitting formulae to the 
surface brightness profiles or as self-consistent dynamical 
models. It is incorrect, however, to infer the full dynam- 
ical properties of the King model solely on the basis of a 
successful fit to the surface density profile. 

The second generically successful set of models was 
empirically derived to reproduce the surface brightness 
profiles of youn g (< 300 Myr) c lusters in the Mag- 
ellanic Clouds ijElson et alJ Il987t lElsonl Il991|) . here- 
after EFF profiles. These models have also been suc- 
cessfully applied to yo ung clusters in external galaxies 
ijWhitmore et alJll999|) . The lack of the developed tidal 
cutoff is associated with the cluster's dynamical youth, 
although profiles of some old Galactic clusters show sim- 
ilar features at low surface bri ghtness levels that ar e 
referred to as extra-tidal stars l|Grillmair et al.l [l995). 
When, or even if, a cluster might evolve from an EFF 
profile to a King profile is unknown and may depend on 
its environment as well as its individual characteristics. 
The lack of a well established tidal limit was confirmed 
for one young LMC cluster (NGC 1866) on the basis of 
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stellar kinematics and more margi nally suggested for tw o 
others (NGC 2164 and NGC 2214: lLupton et all {l989)) 

Whether one type of profile is more appropriate for 
all clusters in the Magellanic Clouds, or whether the ap- 
propriate profile depends on cluster age is also unknown 
because the existing samples of high-quality profiles pref- 
erentially include young clust ers and bec ause they are 
small (18 clusters in the LMC llElsonl 19911) out of an esti- 
mated population of several thousand (B ica et al.1 11999) . 
and none in the SMC) 1 . Previous cluster catalogs have 
identified thousands of clusters in th e Large an d Sma ll 
Magellanic Cloud s (LMC and SMC: IBica, et al l IpO^ : 
iPietrzvhski et all (119981): iBica and Dutra (^OOOD ) over a 
large range of ages llvan den Berghlll981 ). The combi- 
nation of structural parameters and ages is critical in 
unraveling the evolution of these clusters, and by exten- 
sion the processes that drive cluster evolution in general. 
Many issues, s uch as the r elationship between core ra- 
dius and age dElson et al.lll989j) or ellipticity and age 
(jFrenk and Falll 11982) are independent of whether an 
EFF or King profile is most appropriate at large radius. 

Even at the distance of the Magellanic Clouds it can 
often be difficult to resolve individual cluster stars, and 
hence both the quality of the data and the details of 
the model fitting become important. Studies of cluster 
structure fall into two camps. Either they have used 
hi gh quality data for a few clusters, such as the study 
bv iMackev fc Gilmord (|2003albT ) who use Hubble Space 
Telescope images for 10 clusters in the SMC, or they 
have used photograp hic plates, such as the study by 
iKontizas et alJ l)1982|) , which achieve a wider field of view 
but are unable to resolve the inner structure of any clus- 
ter. We pursue a study in the middle ground by adopting 
modern data and algorithms to measure and fit the clus- 
ter luminosity profiles, but use ground based data for 
which determining the structure inside of a few arcsec is 
limited. We present the necessary basic structural pa- 
rameters for 204 SMC clusters. 

This study has two principal purposes. First, we aim 
to determine if clusters, particularly those currently dis- 
solving, are still well described by a King profile, and if 
the older clusters are well described by an EFF profile. 
Second, aside from issues related to the outer profiles, we 
present measurements such as core radius and elliptic- 
ity that are potentially related to cluster evolution and 
independent of the detailed shape of the outer surface 
brightness profile. In addition to the results discussed 
here, these measurements provide a database from which 
future, more detailed studies can select targets and form 
a basis for comparison with other m easurements. For ex- 
ample , cluster ages are presented bv|fiaf elski fc Zaritskvl 
( 2004) and a subsequent study will attempt to place both 
sets of data within a consistent evolutionary model. In 
§2 we describe our cluster catalog, in §3 we present our 
methodology for fitting King and EFF profiles and com- 
pare our results to previous studies, in §4 we develop 
some inferences from the results, and we summarize in 
§5- 

1 Mackcv & Gilmorc 1 2003a b) have recently obtained high qual- 
ity surface brightness profiles of 10 SMC and 53 LMC clusters with 
the Hubble Space Telescope, but they note that their profiles do 
not extend suffici ently far in radius to dist inguish between King 
and EFF models IMackev fc Gilmordl2003al) . 



2. DATA AND CLUSTER CATALOG 

The data come from the Magell anic Clouds Photomet- 
ric Survey (MCPS) described by iZaritskv et all l|1997f) 
and IZaritskv et all l|2002|) . In summary we have ob- 
tained U,B,V, and I drift-scan images of the central 
4.5° by 4° area of the SMC using the the Las Cam- 
pan as Swope (lm) telesco pe and the Great Circle Cam- 
era (jZaritskv et al.l ill 9961)). The stellar p hotometric cat- 
alog is provided bv IZaritskv et all l)2002() and we utilize 
it and the original images in this study. The effective 
exposure time is between 4 and 5 min for the SMC scans 
and the pixel scale is 0.7 arcsec pixel -1 . Typical seeing is 
~ 1.5 arcsec and scans with seeing worse that ^2.5 arcsec 
were rejected. The final catalog contains over 5 million 
stars with at least both B and V measured magnitudes. 
The astrometry for individual stars has a standard de- 
viation of 0.3 arcsec (dervied both from internal errors 
and comparisons wit h a variety of external studies, see 
IZaritskv et"afl l)2002|) '). 

Although extensive cluster catalogs already exist for 
the SMC, we performed our own cluster search to ex- 
amine the completeness of those catalogs over the en- 
tire MCPS survey area. The approach we use is sim- 
ilar to that employe d for a portion of the LMC by 
IZaritskv et alJ (119971) and for th e central ridge of the 
SMC bv IPietrzvhski et~aT] (|1998f) . Specifically we con- 
struct a stellar density image using the stellar photo- 
metric catalog. Within 10 arcsec pixels, we count the 
number of stars with V < 20.5. This particular mag- 
nitude limit is chosen to minimize spurious detections 
and spatial variations due to differing levels of complete- 
ness within and across different scans. The resulting 
stellar density image is then median filtered with a box 
of size 10 such pixels (100 arcsec), and the filtered im- 
age is subtracted from the unfiltered image to remove 
larg e-scale stellar density fluct uations. We use SExtrac- 
tor ijBertin and Arnou"tall996[) to identify peaks in the 
residual image and set the parameters generously so that 
many spurious peaks are found. 

Our need in this particular study is for unambiguous 
clusters for which structural parameters can be mea- 
sured. Most of the effort in building a cluster catalog 
goes into recovering and confirming the most marginal 
objects. As such, some of the clusters published in 
previous studies are quite dubious and various authors 
disagree on the validity of these poor clusters (for ex- 
ample, see the cross-co mparison of the OGLE catalog 
(IPietrzvhski et aj] I1998D with the previous catalog by 
IBica and Dutral (|2000)). Because our aim is to proceed 
with further detailed study of clusters, we focus on defin- 
ing a set of robust clusters. From a cross-comparison of 
our catalog with published catalogs, and then visually in- 
specting the images of candidates, we settle on the list of 
204 clusters presented in Tabled Of the rejected clusters 
from our catalog, many match identifications in other 
catalogs but are either visually unconvincing, too poor to 
offer much chance of further investigation, or embedded 
in strong nebular emission. Therefore, exclusion from our 
list does not necessarily imply that we reject that candi- 
date cluster as a real cluster. The added criteria that we 
have imposed will preferentially exclude young clusters 
within emission regions, those embedded in the denser 
parts of the SMC, and those with a relatively low cen- 



Structural Parameters of SMC Clusters 



3 



Cluster RA DEC Alternate 

Designations 



TABLE 1 
Cluster Designation and Position 



h 24 m 18.67 s 
h 24 m 43.16 s 
h 25 m 26.60 s 
h 27 m 45.17 s 
h 28 m 02.14 s 
h 29 m 55.22 s 
,l 30 m 00.27 s 
h 31 m 01.34 s 
h a9 m 4i n9 s 



-73°59'35.8" 
-73°45'11.7" 
-74°04 / 29.7" 
-72°46'52.5" 
-73°18'13.6" 
-73°41'57.1" 
-73°22 / 40.7" 
-72°20'29.0" 

-79°S4'fin 1" 



K5, L7, ES028-18 
K6, L9, ESO28-20 
K7, Lll, ES028-22 
K8, L12 
HW3 
K9, L13 
HW5 
T.14 




1 0' 

2 0' 

3 0' 

4 0' 

5 0' 

6 0' 

7 0' 

8 0' 

9 0' . -,_ ,« ..... 

10 Q h 32 m 56.26 s -73°06'56.6" NGC152, K10, 

L15, ESQ28-24 



NOTE. — The complete version of this table is in the electronic edition of the 
Journal. The printed edition contains only a sample. 



tral surface brightness. This sample is not complete, nor 
necessarily representative of the overall population, but 
the biases are qualitatively understood and the sample 
is still far larger than any for which detailed structural 
parameter fits are available. 

For the crossidentification presented in Table ^ we 
identify the cl o sest cluster in projection from the 
IBica and Dutral l)2000|) catalog which contains clusters, 
association, and some other sources such as planetary 
nebulae. For matches, the coordinates typically differ by 
a few arcsec, well within the quoted general uncertainty 
of 10 arcsec (|Bica and Dutral l)200(in . 

3. PROFILE FITTING 
3.1. Methodology 

Although we fit both King and EFF profiles to our 
clusters, we place greater emphasis on the King profiles 
for several reasons. First, the full set of King models 
include profiles with extended wings. The largest con- 
centration models we fit resemble a power law profile 
with an effective index similar to the typical EFF pro- 
file (see Figure^. Second, previous studies have only 
demonstrated that the youngest LMC clusters, those 
with age < 300 Myr, are better fit by an EFF profile, 
while over half of our cluster sa mple is estimated to be 
older (Rafclski & Zaritskv 2004). Finally, focusing on a 
single family of models allows for a more straightforward 
comparison across the sample. Our decision to emphasize 
the King profile is certainly open to critique, particularly 
if the data were to favor the EFF profile. Our results 
indicate that this is not the case. If the surface bright- 
ness profiles of young clusters are better fit by softened 
power law profiles, we should find these clusters more 
successfully fit by increasingly concentrated King mod- 
els. However, our analysis shows King models are not 
only statistically satisfactory fits for ~ 90% of the clus- 
ters, but also are generally superior fits than EFF profiles 
for the bulk of the sample (see §3.2). 

The usual approach of fitting a model profile to the 
surface density of a cluster is complicated by the incom- 
pleteness of our stellar catalog that results from crowd- 
ing in the centers of many clusters. We overcome this 
difficulty by fitting models to the surface brightness pro- 
file rather than attempting the difficult and inherently 
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Fig. 1. — A comparison of King and EFF profiles (arbitrary 
units). The short dashed lines represent King models with c = 
0.52, 1, 2.1 and an identical core radius (5), with increasing c going 
from left to right. The solid line represents an EFF profile with 
the same value for a and a value of 7 typical for young clusters, 2.6 
(El son ct al. 1987). The dotted line represents a King model with 
c = 2.1 and a core radius adjusted to result in an outer profile that 
agrees well with that of the EFF model. 



uncertain task of defining and applying an incomplete- 
ness correction. We choose to work with the F-band 
images because they are the deepest of our set, and 
because V is more sensitive to the underlying stellar 
populations than either U or B while being less sensi- 
tive to foreground contaminating stars than /. Further- 
more, u se of V allows us to compare directly to previous 
studies lIKontizas et al.lll982UKontizas fc KontizaslH983t 
iKontizas et al.lll986t iMackev fc Gilmordl2003rjj) without 
concern for possible color gradients. 

Obtaining an accurate determination of each cluster's 
center is the first step in constructing the surface bright- 
ness profile. We use the approximate coordinates given 
in our catalog to extract a 350" x 350" image of the 
cluster. A precise center is determined by calculating 
the luminosity-weighted centroid. To avoid spurious re- 
sults due to contamination by other clusters or bright 
foreground stars, we constrain the calculation to pixels 
within 20" of the approximate center and exclude pixel 
values that are 3000 counts above the mode of the en- 
tire frame (typical mode values are between 100 and 200 
counts). We iterate four times, with each successive run 
using the previous luminosity-weighted center. The de- 
rived centers are stable. For some irregular clusters, we 
adjust the center to match a visually identified center 
that corresponds to central condensation of stars, rather 
than a broader luminosity centroid. This adjustment is 
typically not more than a few arc sec. 

As has been previously noted l|Frenk and Fall ll982). 
many SMC clusters are elliptical. The nature 
of this ellipticity, whether it is fundame ntally con- 
nected to age or luminosity, is debated l)Elsonl 119911 
Ivan den Bergh and Morbevlll984)) . Despite the conclu- 
sion from lElsonl l)1991j) that the asphericity primarily re- 
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suits from unrelaxed stellar concentrations rather than 
any true underlying elliptical light distribution, we opt 
to effectively use elliptical isophotes in determining our 
surface brightness profile. Because we model the sur- 
face brightness of the cluster, subconcentrations should 
be considered both in determining whether a profile is a 
good fit and whether elliptical isophotes are necessary. 

The ellipticity and position angle of each cluster are 
determined by computing the seco nd moments of the 
cluster luminosity distribution (see iBanks et alJ l)1995f) 
for a similar treatment of a small sample of LMC clus- 
ters). We restrict the computation to pixels with fewer 
than 5000 counts above the mode, and to radii between 
3" and 0.7r t , where r t is the fitted tidal radius, or 30" on 
the first iteration of this procedure. These criteria are 
selected with the intention of minimizing the effect of 
background and foreground contamination and ignoring 
the seeing-smoothed inner radii. The ratio of the eigen- 
values of the matrix formed from the second moments of 
the luminosity distribution is the projected ratio of the 
semi-minor, b, to semi-major, a, axis of the cluster, while 
the eigenvector matrix is the rotation matrix of the clus- 
ter. Using these results, we rotate the image so that the 
major axis lies along the horizontal axis and scale the 
image along the vertical axis by a/b. 

We compute the surface brightness profiles using cir- 
cular annuli of 2" width and adopt an uncertainty that 
is given by a of the mean of the pixel values within that 
annulus. We exclude pixels with more than 12000 counts 
above the mean for the annulus to eliminate contamina- 
tion from the brightest stars. We choose different pixel 
thresholds for our calculation of the center, ellipticity, 
and surface brightness because of different sensitivities to 
non-cluster objects. Surface brightness profiles are mea- 
sured out to 200" and include a substantial sky region 
for all clusters. 

The determination of the sky level is key in assessing 
the nature of the outer regions of the surface brightness 
profiles of clusters. As shown in Figure^] an error equiva- 
lent to 10% of the central surface brightness of the cluster 
can convert an intrinsic EFF profile into an apparently 
tidally truncated profile or an intrinsic King profile into 
an apparent one with no tidal radius. Errors in the back- 
ground that are less than 1% of the cluster's central sur- 
face brightness do not qualitatively affect the nature of 
the outer radial profile. It is important to note that this 
is a statistical uncertainty that also affects, but was ig- 
nored by, previous studies. We measure the background 
using an outer annulus that typically extends over the 
last quarter of the radial range of the extracted images 
(r > 150"), although we adjust the boundaries of the sky 
region interactively for large clusters or for those with 
contaminat ing stars. Because youn g stars are spatially 
correlated {Harris fc Zaritskvl llQQQI. going to larger and 
larger radius to measure the background is not an ap- 
propriate solution, particularly for the younger clusters 
where the spatial correlations are the strongest. More- 
over, there is no reason to expect the background at large 
radii from the cluster to be the same as the background 
in which the cluster resides, especially if there are many 
young stars in the region. Unfortunately, adopting back- 
ground annuli near the clusters can often result in back- 
ground measurements that are affected by local fluctua- 
tions. The determination of the background is therefore 



-J 











i i „.,!,,,, i \ \ 


L+10% \ 




L-10% J 




j 1 1 1 1 1 1 1 


+t+f 


:l 1 1 1 1 1 1 1 1 1 1 1 1: 




L+i% 

= 


ilF 


E i f 





jj 1 1 1 1 | 1 u 










L+10% 1 




L-10% J 




:l 1 1 1 1 1 1 


\mi 


:l 1 1 1 1 1 1 1 1 1 1 1 1; 




L+1% 




L-US^^y _ 



log 10 (Radius) 

Fig. 2. — A comparison of the effect of background subtrac- 
tion errors on both representative King and EFF profiles. Upper 
set of panels illustrate the effect of oversubtracting (+) and un- 
dersubtracting (— ) the background by a percentage of the central 
brightness of the cluster on EFF profiles (the original profile is rep- 
resented by the solid line in the upper left panel). The dashed line 
represents the comparison King profile in the upper panel and the 
comparison EFF profile in the lower panel. The lower panels are 
analogous except for illustrating the effect of background errors on 
a King profile. 



an inherently poorly solved problem. We opt against in- 
cluding the background as a free parameter in the model 
fitting described below because we find that including 
radial bins in which the background dominates allowed 
us to improve the statistical measure of the fit, although 
it often introduces systematic errors in the fitted clus- 
ter profile. This problem arises because the photometric 
properties of the area around the cluster are not neces- 
sarily identical to those underlying the cluster, in partic- 
ular the background is not uniform. A statistical fitting 
method attempts to fit the background at the expense 
of fitting the cluster and results in fits that are visu- 
ally inferior in many cases. Including only radial bins in 
which the cluster dominates results in fitting degenera- 
cies due to the lack of a strong direct constraint on the 
background level. Instead, we determine the background 
independently from a sky annulus at moderate radius, 
which is similar to the procedure adopted by EFF. The 
inner boundary of our sky region is typically about 8 
times larger than the radius that includes 90% of the 
cluster light. Our fits to the background, and the com- 
parison benchmarks of 1 and 10% errors, are shown in 
Figure|s3 The Figure illustrates that there are some clus- 
ters for which the scatter in the background exceeds the 
10% benchmark and hence sky uncertainties play a sig- 
nificant role, while there are others for which the scatter 
is less than the 1% benchmark and hence sky uncertain- 
ties are sufficiently small that the outer cluster profile is 
well measured. 

We fit King and EFF profiles to each of our clusters 
via chi-squared minimization. The standard King model 
is des cribed by a combination of three parameters {Kind 
1966). We choose to define each profile by its central 
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Fig. 3. — Background fits for the first 25 SMC clusters (numbers in lower right are the cluster designation). The range of horizontal axis 
is from to 200 arcsec for all clusters. The vertical scale is unique for each cluster and corresponds to ± 10% of the selected background 
value. The sky used in our fitting is shown by the solid line. The short dashed line mark the range of error corresponding to 1% of the 
central surface brightness of the cluster and the long dashed lines correspond to 10%. We show a representative sample here and present 
all the clusters in the electronic version of the Figure. 



surface brightness, So, core radius, r c , and concentra- 
tion, c. Other common parameters, such as tidal radius, 
can be derived from this set (ie. c = logio(r t /r c )). In 
addition to these three parameters, our profiles include 
the background surface brightness discussed above. We 
explore a parameter space that ranges in Eo from 0.25 
to 1.75 times the mean E measured within the central 
4" of the cluster in increments of 0.025, in r c from 1" to 
101" in increments of 0.5", and in c from 0.001 to 2.1 in 
increments of 0.01. The concentration range corresponds 
to the range over which King's original e mpirical models 
match his theoretical models as found bv lWivanto et all 
(1985). The fitting is iterated once. Error bounds on 
each best fit parameter corresponds to the range of that 
parameter among models that cannot be rejected with 



more than 67% confid ence. 

As mentioned above . lElson et aTl l|1987j) noted that the 
surface brightness profile, /Lt(r), of young clusters appear 
to be better described by a core plus power law profile, 
fi(r) oc (1 + r 2 /a 2 ) -7 / 2 , where a and 7 are free parame- 
ters. We fit the EFF models our data in the same manner 
as we fit the King profiles, over the same radial range, 
by varying the central density (Eo), the scale length (a), 
and the power-law index (7). Note that a is the scale 
length, rather than the core radius, which can be cal- 
culated using r c = av / 2 2 / 7 — 1. We explore a parameter 
space that ranges in Eo from 0.25 to 1.75 times the mean 
E measured within the central 4"of the cluster in incre- 
ments of 0.025, in a in from 1" to 101 " in increments 
of 0.5", and in 7 from 1.5 to 7.0 in increments of 0.05. 
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Fig. 4. — V-band surface brightness profiles of the first 25 SMC clusters and profile fits. These profiles include the background flux. The 
range of horizontal axis is from to 200 arcsec for all clusters. The vertical scale is unique for each cluster. The central surface brightness 
of each cluster fit is given in Table |21 The best fit King profile is shown by the solid line. The numeral in the upper right of each panel 
provides the cluster identification. We show a representative sample here and present all the clusters in the electronic version of the Figure. 



We present the best fit King profiles in Figure 0] and 
both King and EFF profiles in Figure [SJoverlayed on the 
data. The key difference between the two Figures is that 
Figure 0] shows the un-background subtracted data so 
that all of the data could be plotted on the logarithmic 
scale. Note however, the profiles plotted in this man- 
ner do not resemble the standard King models, although 
they are. In Figure we again plot the background- 
subtracted profiles, but here we stop plotting data and 
models once a radial bin has a negative background sub- 
tracted value because the logarithm is undefined. Plot- 
ting only data greater than the background, i.e. those 
with positive background subtracted values, results in a 
misleading profile. The parameter values for the fits are 
presented in Tabled The la confidence intervals are de- 
fined using x 2 , and are therefore undefined in the cases 



where even the best fit model has \ 2 greater than the 
corresponding value for the available degrees of freedom. 

3.2. Results 

Of the 204 clusters in our catalog, only 44 are poorly 
described by a King profile with greater than 67% confi- 
dence. Given the size of our sample however, we expect 
from statistical considerations that a substantial number 
of clusters that are truly well described by a King pro- 
file will be rejected at this level. We calculate that only 
19 clusters are fit sufficiently poorly that for each indi- 
vidual cluster the King profile can be rejected as an ac- 
curate description of their surface brightness profile (this 
limit corresponds to rejecting the profile with 99.5% con- 
fidence) . Because both the King and EFF profiles can fit 
clusters with extended radial surface brightness profiles, 
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0.5 1 1.5 2 0.5 1 1,5 2 0.5 1 1.5 2 0.5 1 1.5 2 0.5 1 1.5 2 

log 10 (Radius) [arcsec] 

Fig. 5. — The V-band cluster profiles, after background subtraction, are compared to both the best fitting King (solid line) and EFF 
profiles (dashed line). The radial plots are truncated at the innermost radius for which the counts minus the background are less than zero 
to avoid only plotting values for which the logarithm exists. This omission results in some profiles that appear to be systematically poor fits 
because the fit includes data that are omitted from the Figure. Error bars have been omitted for clarity. All panels are identically scaled 
and the numeral in the upper right corner provides the cluster identification. All cluster profiles are available in the electronic version of 
this Figure. 



the profile itself cannot be used to infer extratidal stars 
(i.e. stars beyond the cluster's Roche lobe). While the 
extended wings in the EFF profile are generally inter- 
preted to arise from that type of mass loss, the extended 
wings in a high-concentration King model are bound to 
the system. Determining which physical interpretation 
is correct will require either direct dynamical measure- 
ments, such as a measurement of the cluster mass and the 
velocity dispersion of stars, or indirect evidence, such as 
a lack of old clusters with such wings, which would imply 
that there has been dynamical evolution of such clusters. 
In this study, we reach conclusions on the relative mer- 
its of the fitting functions, not on the dynamical state 
of the clusters. As we have illustrated before (Figure 
describing a surface brightness with a King profile does 



not necessarily imply the presence of a dramatic tidal 
edge. 

We find that our determination of the tidal radius is 
highly uncertain. This uncertainty arises from the strong 
coupling of the estimated r t to the determination of the 
background. Similarly, the uncertainty on the concentra- 
tion, because it depends on r t , is too large to be meaning- 
ful. Because we desire a robust measure of cluster size, we 
pursue an alternative definition of size determined by the 
enclosure of a significant fraction of the total cluster light, 
rather than all of the light. We choose to tabulate the 
radius inside of which 90% of the total cluster luminosity 
of the fitted profile is contained, rgo. Even though rgo, 
by definition, still contains the bulk of the cluster light, 
the associated uncertainty is typically half that of r t . We 
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Fig. 6. — V-band images of the first 25 SMC clusters. The 
angular size of each image corresponds to the central quarter of the 
area for which the profiles in Figure |1] are de rived. Images of all 
the cl usters in the four filters are available in Rafclski & Zaritsky 
(2001). 
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Fig. 8. — Comparison of our measured core and tidal radii to 
previously published values. Our core radii are compared to those 
presented in Kontizas et al (1982,1983,1986) in the upper left panel 
and those measured by Mackey & Gilmore (2003b) in the upper 
right. The core radii of Kontizas et al and Mackey & Gilmore 
are compared at bottom left. The bottom right panel features a 
comparison on our tidal radii to those published by Kontizas et al. 
Clusters at (135,26) in the upper left and (55,207) in the bottom 
right were excluded for clarity. 
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Fig. 7. — Comparison of rgo values obtained with the King and 
EFF models. The solid line is the 1:1 line. 



avoid going to even smaller radii because then connection 
to the cluster size becomes increasingly tenuous and de- 
pendent on the cluster profile. Because of the ambiguity 
of which model profile to use, one may wonder whether 
defining rgo using one profile or the other leads to large 
differences in the calculated value of rgo . We show in Fig- 
urcQa comparison of the values of rgo obtained using the 
best fit King or EFF model. We find that although there 
is some systematic difference between the two that can 
be as large as 25%, the rankings of clusters according to 
size (and hence any correlation one might examine using 



rgo) is maintained to large degree regardless of the model 
one adopts. The best fit parameters for each candidate, 
along with the la uncertainties, are presented in Table 
The listed magnitude is calculated by integrating the 
fitted King profile. The complete set of profiles, images, 
and structural parameters is presented in the online ver- 
sion of these Figures and Tables. 

Our measurements of structural parameters constitute 
the largest such survey of SMC cluster profiles, so the 
availability of comparison data is limited. There ex- 
ist only four published studies of SMC cluster profiles, 
three of which are by a single group of authors and 
should be considered a sing le survey. 
1 1982 1) . iKontizas & Kontizas! IpS.U and 



Kontizas et al 



Kont izas et al 



11986), hereafter collectively referred to as Kontizas et 
al., use photographic plates from the U.K. Schmidt and 
Anglo-Australian tel escopes in Australia to me asure 67 
cluster profiles, while lMackev fc Gilmord l)2003b|) present 
10 profiles measured from images obtained with the 
Hubble Space Tele s cope. All 10 clusters observed by 
Mackev & Gilmore (2003b) appear in Kontizas et al., 
but most were observed using different filter passbands. 
Although we do not know whether the structural pa- 
rameters are strongly dependent on color, a conservative 
comparison should be limited to parameters determined 
from images taken with similar filters. An examination 
of Kontizas et al.'s observations in a variety of filters 
and emulsions do suggest that a cluster's profile may de- 
pend on the wavelength in which it is observed, although 
this issue is not specifically investigated. We restrict our 
comparison to the smaller set of 28 clusters observed in 
a filter and emulsion combination closely corresponding 
to our l/-band observations. Our cluster sample includes 
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TABLE 2 

Structural Parameters' 1 



Cluster Core Radius [pc] 90% Light Radius [pc] My So,K e xl K EFF Scalelength b [pc] 7 S , E 

low best high low best high [cts/ low best high low best high [cts/] 

■3 1 dT , /- > cci/-"2 1 



1 


0.61 


0.61 


0.92 


11.58 


15.56 


17.34 


18.7 


60 


0.30 


0.99 




0.76 






6.30 




421 


1. 


2 


2.44 


3.97 


7.33 


12.17 


18.06 


20.84 


12.9 


595 


0.17 


0.77 


1.07 


3.35 


8.37 


1.50 


2.20 


4.65 


659 


0. 


3 




2.75 






11.30 




14.5 


314 


0.02 


1.25 




2.59 






2.45 




335 


V 


1 


5.19 


10.69 


18.94 


11.54 


12.62 


16.69 


13.8 


133 


0.18 


0.81 




13.01 






6.95 




111 


1. 


5 


1.53 


3.97 


12.52 


6.43 


11.90 


19.20 


14.8 


189 


0.28 


0.68 


1.22 


8.52 


10.96 


1.50 


6.90 


7.00 


167 


0. 


6 


0.61 


3.36 


14.97 


10.39 


18.31 


24.30 


16.9 


145 


0.37 


0.79 


0.61 


1.37 


12.94 


1.50 


2.75 


7.00 


135 


0. 


7 


9.16 


9.16 


9.77 


9.57 


9.57 


9.62 


14.5 


139 


0.56 


1.10 


3.51 


6.55 


11.87 


1.85 


3.40 


7.00 


153 


1. 


8 


0.92 


0.92 


1.53 


7.08 


11.78 


14.98 


14.7 


1259 


0.38 


0.87 


0.61 


0.76 


1.83 


1.80 


2.20 


3.10 


1835 


0. 


9 


2.75 


4.89 


12.52 


6.99 


11.24 


12.81 


14.8 


116 


0.37 


0.70 


3.20 


8.68 


10.50 


2.25 


7.00 


7.00 


116 


0. 


10 


6.72 


9.16 


14.36 


19.82 


22.40 


24.11 


12.1 


474 


0.23 


0.79 




18.57 






6.30 




399 


1. 



NOTE. — The complete version of this Table is in the electronic edition of the Journal. The 

a The ltr confidence intervals arc defined using x 2 , and arc therefore undefined in the cases 
of freedom. 

^ The core radius, when defined as the radius for which the surface brightness drops to half 



printed edition contains only a sample. 

where even the best fit model has x 2 greater than the corresponding value for the available 
its central value, can be calculated using t c = oV 2 2 /7 - 1. 



25 from this restricted set. 

Our measurements of r c a nd r± are compared t o those 
made by Kontizas et al. and lMackev fc Gilmord ((2003b) 
in Figure[5| Our core radii are systematically larger than 
those measured by Kontizas et al. This result is con- 
sistent with Kontizas et al.'s prediction that their core 
radii are underestimated, a conclusion also reached by 
iMackev fc Cilmorel l)2003btL Given that Kontizas et al. 
could not resolve the core in any of their clusters, it is 
not surprising that there is di sagreement with both ou r 
measurements and those of IMackev fc Gilmord (J2003b). 
More illuminat ing is the comparison of rv measurements 
with those of IMackev fc Gilmord l)200 3b ). who are able 
to resolve the inner regions of the clusters. Here we find 
that our results are statistically consistent, although with 
some slight bias to overestimate r c by ~ 30%. We sus- 
pect that our core radii are systematically larger because 
we are unable to resolve the inner cluster structure to the 
degree possible in the HST observations. 

Comparison of tidal radii is more uncertain both be- 
cause there is a smaller set of data with which to 
compare to and the lack of modern techniques used 
in constructing sm a ller co mparison set. Although the 
IMackev fc Gilmorel l)2003bf) observations do not extend 
to large enough radii, Kontizas et al's do. Figure |H1 re- 
veals that our tidal measurements are consistent with, 
although ~ 20% systematically smaller than, those re- 
ported by Kontizas et al. Given the lack of precision with 
which tidal radius can be fit (see H3.2(l . this difference is 
not surprising, especially given the significant difference 
in the quality of the data and analysis techniques between 
our study and that by Kontizas et al. We conclude that 
our model fits are in rough accord with previous fits, but 
given the size of the scatter seen in Figure [S] and our cal- 
culated uncertainties for r t we do not include r t values 
in Table □ 

We now return to the issue of the appropriateness of 
the King model. There are two fundamental issues: 1) 
does the King profile provide an acceptable description 
of the surface brightness profile and 2) does the King 
model provide an acceptable description of the stellar 




King Model x Z 

Fig. 9. — The xl values for the best fit King and EFF profiles. 
The line is the 1:1 line. Only clusters for which the fits have xl < 3 
are shown. Those clusters with background values determined to 
better than 1% of the central cluster surface brightness are plotted 
with larger open circles surrounding the central filled circle. 



dynamics. With our data, we are only in a position to 
address the former of these. Indeed, dynamical consider- 
ations would suggest, as noted above, that King models 
are unlikely to be the correct dynamical description of 
the younger clusters. For the bulk of the clusters there 
are only subtle differences between the King and EFF 
profile that mostly occur near the outer edge of the clus- 
ters, in the region where the background uncertainties 
may dominate. The fit in this region is not shown in 
Figure because we stop plotting once the background 
subtracted flux value is less than zero, and is hence unde- 
fined in the units plotted. However a comparison of the 
\t values, Figure El demonstrates that while for any in- 
dividual cluster the is typically not large enough to 
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Fig. 10.— The A\l values vs. age. We include only clusters for 
which either the King or EFF models fit with xt <1 an< i < 2. 



rule out one model relative to the other, the King model 
systematically fits these cluster profiles better. To ad- 
dress the importance of the background subtraction to 
this conclusion, we have highlighted in the Figure the 
clusters for which the background is determined to bet- 
ter than 1% of the cluster's central surface brightness (an 
uncertainty level below which the cluster profiles do not 
appear to be grossly disturbed by incorrect background 
subtraction, see Figure EJ) . These high quality surface 
brightness profiles follow the mean trend. In the cases 
where the King model does significantly better than the 
EFF models and the cluster is a rich cluster with a well 
determined sky brightness and overall profile, the differ- 
ence in the models is clearly in the outermost regions 
(see for examples clusters 85, 86, 96, 143, 174). In other 
cases, the cause for the preference of the King model is 
more subtle and ambiguous. It is presumably due to a 
cumulative effect in the fit over the range between core 
and tidal radius. 

To examine whether one particular model is favored 
for clusters of a particular age, for example EFF mod- 
els for the youngest clusters, we have plotted the dif- 
ference in y?, value s as a function of age (ages from 
iRafelski fc Zaritskvl l|2004|) 'l in Figure [TBI for fits where 
either the King or EFF model had xt < 1- There is 
no trend evident that favors one particular model at a 
specific age. The general preference for King models 
(Axl > 0) is seen at all ages. Of course this result 
doesn't argue against the effectiveness of the EFF pro- 
file for some clusters, nor against the existence of some 
clusters that are not tidally limited. Nevertheless, the 
situation is clearly more complicated than an expecta- 
tion that all young clusters should have power-law sur- 
face brightness profiles. 

4. DISCUSSION 

4.1. Cluster Subpopulations 

Even those clusters that are well-fit by a King profile 
are not drawn from a single parameter population. Al- 



FlG. 11. — The distribution of cluster concentrations in the SMC. 
Note that rgo as opposed to tidal radius is used to calculate the 
concentrations for this Figure, while the standard definition of con- 
centration is fit in the model. The axes are set so that the entire 
range of allowable values of concentration is shown. 



though the core radius, r c , and rgo are correlated at the 
92.7% confidence level (Spearman rank test; see Table 
EJ) , clusters distribute themselves over the entire range of 
rgo /r c that we allow in our fit ting space given the th eoret- 
ical constraints described bv lWivanto et al.l l|1985|) (Fig- 
ure lll|) . This implies that a variety of factors influence 
the observed structural parameters. Possible influences 
include tidal effects, age (the internal dynamical evolu- 
tion of the cluster), and initial conditions (substructure, 
velocity anisotropics, initial stellar mass function, angu- 
lar momentum, and chemical abundance). These are no- 
toriously difficult to disentangle, but the large, uniform 
sample presented here offers a new opportunity to search 
for clues. 

We find King profiles to be a statistically acceptable 
fit to clusters that range over a factor of over 40 in r c , 
10 in rgo, and 10 3 in So- This finding is somewhat sur- 
prising in light of the fact that the vast majority of clus- 
ters are inferred to be in the process of dissolution or 
disruption on the basis of the cluster age distribution 
ijBoutloukos fc Lamersll2003t IRafelski fc Zaritskvl |200|) 
and many have had insufficient time to dynamically re- 
lax. Therefore, the fits to some clusters are almost cer- 
tainly not physically meaningful, but rather are statisti- 
cally allowed given the large observational uncertainties. 
Nevertheless, even in these cases our measurements do 
reflect the size of the cluster and its core. Although the 
large majority of SMC clusters are adequately fit by sim- 
ple King profiles, and hence demonstrate that measuring 
surface brightness profiles alone is not a promising av- 
enue to identifying dissolving clusters, there are several 
populations of clusters that merit further discussion and 
may be showing signs of their dynamical evolution. 

4.1.1. The Ring Clusters 

A subset of our cluster luminosity profiles contain a 
significant bump in the luminosity profile relative to the 
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Fig. 12. — The distribution of cluster ellipticities in the SMC and 
MW. The SMC distribution is depicted in top panel while that of 
the MW appears in the bottom panel. The MW data come from ? 
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Fig. 13. — V-band images of the 25 most elliptical clusters in the 
sample 



fitted King or EFF profile. Of those, visual inspection 
suggests that about 12 out of the 204 have a bump that 
occurs beyond the central 8 arcsec and is not apparently 
due to a few luminous stars (clusters 12, 19, 39, 102, 
115, 148, 164, 179, 183, 186, 196, and 203 in Figure 
. Although the apparent ring features in some of these 
clusters may be due to statistical fluctuations, there are 
at least some where the central deficit of stars is quite 
evident (see clusters 19, 39, 148, and 186 in Figure HTjl 
and the xt values, for these clusters, are larger than ex- 
pected for random fluctuations in a sample of this size 
with 99.5% confidence for either the King or EFF pro- 
files (see §3.2). We discount similar features at small ra- 
dius, for example in clusters 3 and 4, because they could 
arise from slight centering offsets or random fluctuations 
within the small inner annuli. We also have not included 
some clusters in which the profile bump appears to be 
due to bright contaminating stars (such as cluster 82). 
Nevertheless, there may be many more clusters in the 
sample than the 12 listed in which the bump is a real 
feature. One of the best examples of this phenomenon is 
cluster 19 (Figures and 117(1 . This feature may be a 
signal of a cluster's dissolution because such a configu- 
ration appears to be difficult to maintain in equilibrium. 
Further dynamical modeling is required to confirm this 
supposition. Aside from dynamical issues, it is also the 
case that because we fit to the luminosity profile rather 
than the stellar density profile, we are susceptible to vari- 
ations in the mass (or number) to light ratio. The degree 
to which the underlying light profile departs from a King 
or EFF profile may not indicate a similarly strong devi- 
ation in the mass profile. 

4.1.2. The Ellliptical Clusters 

The distribution of cluster ellipticities is presented in 
Figure To explore the origin of the apparent broad 
distribution relative to that of the Ga lactic clusters (data 
from the compilation of lHarrisI l|1996|^ . we search for sig- 



nificant correlations between ellipticity and other cluster 
properties. For example, elli pticity has been fo und to 
correlate inversely with age l(Frenk and Fall 1982(1 , which 
suggested either evolutionary changes in the structure of 
clusters or changes in the conditions of formatio n. In 
a diffe rent interpretation of the same correlation. lElsonl 
(1991) suggested that subclumps at various stages of 
merging may produce the impression of large ellipticity in 
LMC clusters. This conclusion was based on the inverse 
correlation between age and ellipticity and the interpre- 
tation of various features in clust er profiles. However, 
Ivan den Bergh and Morbevl l(1984(l find that although el- 
lipticity and age are indeed inversely correlated, that re- 
sult is a byproduct of the inverse correlation betw een el- 
liptici ty and luminosity, a claim not addressed bv lElsonl 
l(1991() . However, it is quite difficult to work in terms of 
luminosity since it depends both on mass and age, both 
of which might help determine the ellipticity. 

In Figures El an d El we present the images of the 
25 most and least elliptical clusters. The most ellipti- 
cal clusters are in general also highly irregular, as well 
as of low luminosity. The least elliptical clusters are in 
general, although not exclusively, richer systems. While 
there are comparably poor, low luminosity systems that 
are both of high and low ellipticity, the luminous systems 
appear to uniformly avoid the upper tail of the ellipticity 
distribution. Two examples of rich clusters with moder- 
ately high ellipticities (~ 0.2) are presented in FigureEl 
These are among the most elliptical of the rich clusters 
in the SMC and interestingly correspond to the highest 
ellipticities measured for Galactic cl usters. 

We use the age dating results of IR afelski fc Zaritskvl 
(2004) to examine whether age or mass is the primary 
driver of ellipticity. We hav e converted the age and lu- 
minosity measurements from IRafelski fc Zaritskvl ( 2004) 
for the Z = 0.004 m odels into a mass by u sing the same 
Starburst99 models ((Leitherer et alJll999(l used to esti- 
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Fig. 14. — V-band images of the 25 least elliptical clusters in the 
sample 



mate the ages. In Figure [TBI we plot ellipticity vs. mass 
for different age bins. The three bins are defined such 
that the youngest bin contains clusters whose 90% upper 
confidence age is less than 500 Myr, the middle bin con- 
tains clusters whose best fit age is between 1 and 3 Gyr, 
and the oldest bins contains clusters whose 90% lower 
confidence age is greater than 3 Gyr. The mean elliptic- 
ity of clusters in the three age bins are 0.31, 0.27, and 
0.32 (from youngest to oldest), and so we find no corre- 
lation between age and ellipticity. On the other hand, all 
three age bins show some sign of a correlation with mass. 
All the data together results in a significant detection of 
an inverse correlation (Spearman rank correlation con- 
fidence levels of 99.1%). Results for the individual age 
bins suggest similar inverse correlations, but the correla- 
tion is only significant in the intermediate age bin (97.5% 
confidence). Examining the joint panel illustrates that it 
is not possible to form large samples of clusters of similar 
mass with which to examine trends with age. For exam- 
ple, all but seven of the young clusters have masses less 
than 5000 M©, while none of the older clusters have such 
low masses. Selecting clusters by luminosity mixes lower 
mass young clusters with higher mass old clusters and 
so results in a comparison of more massive old systems, 
with less massive young systems. While our results show 
that the dependence of ellipticity on mass is more easily 
identified than any with age, there is a relative deficit of 
young clusters with very low ellipticity, suggesting that 
age, and dynamical evolution, may still play a role. 

4.1.3. The Irregular Clusters 

In FigureH7lwe show the images of the 20 clusters with 
the highest xt values relative to the King model fitting. 
These also have high xt values relative to the EFF model 
fitting. Several of these clusters are the "ring"-type clus- 
ters described above (19, 39, 148, 186), which explains 
why these are poorly fit by the either profile. Others 




Fig. 15. — V-band images of cluster 98 and 197. 
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Fig. 16. — Cluster ellipticity vs. mass for three different age bins. 
The left panel shows all the data, where the youngest clusters at 
represented with circles, the intermediate age clusters with squares, 
and the oldest clusters with triangles (see text for defining ages). 
The right panels separate the three age populations for clarity and 
rescale to maximize the plot for the relevant population. 



have very high central luminosities (70, 71, 181). The 
remainder appear to suffer from luminosity fluctuations 
caused by a few stars, which are perhaps unassociated 
with the cluster. Our principal conclusion is that with 
the exception of the ring clusters, we find no system- 
atic failure mode for the King and EFF profile fitting. 
Even in cases where xt ^ s iar S e J the fits appear to be 
sufficiently good that one can ascribe the mathematical 
failure to such systematics as contamination. 

4.2. Correlations and Comparisons 

We present a set of Spearman rank correlation mea- 
sures among the cluster parameters we considered (Table 
|3J). Such an exercise can often lead to spurious correla- 
tions that are in fact driven by hidden variables, but it 
can also provide clues to fundamental physical drivers. 
We discuss a few of the statistically significant corre- 
lations that we find interesting. First, we find an in- 
verse correlation between core radius and central surface 
brightness. Such a correlation would appear to arise nat- 
urally among a set of clusters that are dissolving. Second, 
we find a correlation between the central surface bright- 
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Fig. 17. — V-band images of the twenty clusters with the largest 
X^'s. From upper left across and then down these are clusters 71, 
89, 19, 17, 61, 155, 88, 186, 39, 15, 70, 145, 87, 181, 114, 53, 148, 
81, 77, and 143. 
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ness and the background level. This could arise from 
the practical difficulty in identifying low-surface bright- 
ness clusters in denser environments, but could also be 
due to more effic ient disruption of such clu sters in higher 
density regions (Gnedin & Ostriker 1997). These cor- 
relations are all internally consistent with a model that 
has clusters preferentially disrupting in higher density re- 
gions. One can only proceed so far on the basis of these 
correlations. Two key inputs are necessary to continue: 
1) ages to disentangle evolutionary effects and 2) models 
that make quantitative predictions for these correlations. 
We are c urrently working on Joinin g the available age 
estimates (iRafelski fc Zaritskvll2004) with model predic- 
tions to address these questions. 

To complement the internal comparison approach, we 
also compare the cluster properties to those of Milky Way 
cl usters (data fo r Galactic clusters from the compilation 
of Harris (1996)). We present a comparison of the core 
radii of Milky Way and SMC clusters in Figure CHI The 
distribution of core radius size is much broader in the 
SMC than in the MW. This difference is either due to the 
preferred dissolution of young, low-concentration clusters 
(both in the SMC and MW) and the preponderance of 
young clusters in the SMC, or to the preferred forma- 
tion and survival of lower-concentration clusters in the 
SMC. The distribution of core radii for clusters with age 
> 1 Gyr (Figure ITH|l shows that the older population is 
just as broad in its distribution of core radii and demon- 
strates that long-lasting, low-concentration clusters are 
more prevalent in the SMC than in the MW. Similar re- 
sults were found for a smaller SMC cluster sample, as 
well as cluste r samples for the LMC, F ornax, and Sagit- 
tarius dSphs i|Mackev fc Gilmordl2T)0l . 

5. SUMMARY 

We present structural parameters for 204 stellar clus- 
ters in the Small Magellanic Cloud derived from fitting 
King and Elson, Fall, and Freeman (EFF) profiles to 
the V-band surface brightness profiles as measured from 
the Magellanic Clouds Photometric Survey images. We 
find that King profiles do a remarkable job at fitting the 
vast majority of the clusters, which is particularly sur- 
prising given that 99% of the clusters will dissolve over 



Fig. 18. — The distribution of core radii in the SMC and MW. 
The SMC distribution is depicted in top panel while that of the 
MW appears in the bottom panel. In the top panel we also show 
the distribution for SMC clusters with age > 1 Gyr in the shaded 
historgram. The MW data come from ? 



the n ext few hundred million years ijRafelski fc Zaritskvl 
2004) and that many have had insufficient time to dy- 
namically relax. We find that the King profiles system- 
atically fit slightly better than the EFF model for most 
clusters in our sample, independent of cluster age, but 
that on a cluster by cluster basis it is difficult to rule out 
one type of model or the other. The only systematic de- 
viation in the surface brightness profiles that we identify 
is in a set of clusters that lack a central concentration, 
and which we have named as "ring" clusters. In agree- 
ment with previous studies, we find that the clusters in 
the SMC are significantly more elliptical than those in 
the Milky Way. We are able to identify an inverse cor- 
relation between cluster mass and ellipticity, but find no 
evidence for one between age and ellipticity. The search 
for the latter is hampered by the lack of a sizeable sam- 
ple of clusters with similar masses and a range of ages. 
We identify several correlations (central surface bright- 
ness vs. local background density, size vs. distance) that 
an be used to constrain models of cluster evolution Much 
of the interpretation of the empirical results we present 
rests on having measured ages and a detailed model for 
the formation and rapid evolut ion of the c l uster popula- 
tion (extending the models of iBekki et all l|2004f) \ This 
current sample of structural parameters provides the ba- 
sis for detailed studies of the evolution of the cluster sys- 
tem of the Small Magellanic Cloud. 
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TABLE 3 

Spearman Rank Correlation Probabilities 



Parameter r c rgo ego Central Den. Luminosity Ellipticity Background 



r c 




.927 


-1 


-1 


1 


.741 


-.957 


rgo 


.927 




1 


-.434 


1 


-.678 


.369 


C90 


-1 


1 




1 


-.613 


-.659 


.951 


Central Den. 


-1 


-.434 


1 




1 


-.974 


.998 


Luminosity 


1 


1 


-.613 


1 




-.991 


.909 


Ellipticity 


.741 


-.678 


-.659 


-.974 


-.991 




-.874 


Background 


-.957 


.369 


.951 


.998 


.909 


-.874 





NOTE. — Probability of correlation between listed cluster properties computed with Spearman rank analysis. Negative 
sign denotes anti-correlation. 
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